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We present a technique for the successive restoration of the branches of the one-dimensional 
strip brightness distribution across a quasar's accretion disk via the analysis of observations 
of high magnification events in measured fluxes from the multiple quasar images produced 
by a gravitational lens. Hypothesizing these events to by associated with microlensing by a 
fold caustic, the branches of brightness distribution are searched for on compact sets of non- 
negative, monotonically non-increasing, convex downward functions. The results of numerical 
simulations show that the solution obtained is stable against random noise. Analysis of the light 
curves of high magnification events in the fluxes from components C and A of the gravitational 
lens QSO 2237+0305, observed by the OGLE and GLITP groups, has yielded the forms of the 
strip brightness distributions across the accretion disk of the lensed quasar. The resulting sizes 
of the accretion disk are in agreement with results obtained earlier via model-fitting. The form 
of the brightness distribution is consistent with the expected appearance of an accretion disk 
rotating around supermassive black hole. 



1. Introduction 



According to our current understanding, the main source of the energy of quasars and other 
active galactic nuclei (AGN) is the accretion onto a supermassive black hole. The properties of 
the accretion disk that is formed determine to an appreciable extent the observed characteristics 
of these objects [1]. One way to study accretion processes is via spectral analyses, including 
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investigations of the X-ray lines profiles and the reverberation mapping or echo tomography of 
the disks [2]. Such studies have yielded most of the currently available information. 

Another independent way to obtain information about the accretion processes is to investigate 
the spatial structure of the disk, which requires observations with very high angular resolution, 
exceeding a microsecond of arc. In spite of the seemingly unrealistic smallness of this value, a 
proposed space X-ray interferometer project [3] would, in fact, provide this resolution. However, 
a similar resolution can be obtained in the visible using observations of gravitational lens 
systems: intervening galaxies that give rise to multiple images of distant quasars. Microlesing 
by stars in the lens galaxy produces a random field of caustics, which can lead to a high 
magnification event (HME) in the measured fluxes from the images when a caustic crosses the 
accretion disk of the quasar [4,5]. The most probable type of microlensing occurs when a fold 
caustic intersects the disk, in which case the observed light curve contains information about 
the one-dimensional strip brightness distribution across the accretion disk in the direction of 
the local normal to the caustic. Restoration of this strip distribution from observational data 
involves solution an ill-posed inverse problem, and requires the use of special algorithms that 
are stable against random noise. This problem was first considered by Griegar et al. [6], who 
used the Tikhonov's regularization method. 

The more interesting problem of restoration of the radial brightness distribution in the locally 
co-moving frame of the accretion disk was first analyzed by Agol and Krolik [7], who also used 
the Tikhonov's regularization method. This problem is significantly more difficult, and requires, 
in addition to the introduction of a number of free parameters to describe the microlensing 
geometry, the consideration of relativistic effects. For observations obtained in a single photo- 
metric band, it is not possible to take relativistic effects into account without a disk model that 
describes how the radiation intensity of a disk area element depends on direction and frequency. 
The attempt of Mineshige and Yonehara [8] to restore the radial brightness distribution in the 
accretion disk seen by an external distant observer was not successful in this sense, since their 
solution assumed circular symmetry, which will not be the case in the presence of relativistic 
effects. In spite of the possibility of deriving the radial brightness distribution in an accretion 
disk in a locally co-moving frame, the restoration of the strip brightness distribution in a frame 
of a distant observer remains an important problem. This is true because the form of the strip 
brightness distribution can be restored using minimum a priori information without knowledge 
of the caustic parameters, distances to the source and lens, or the relative velocity of the source 
and lens. There is no need to introduce any model constraints on the properties of the lensed 
source. 

The gravitational lens QSO 2237+0305, which is also called Huchra's lens or the Einstein cross, 
is the best know representative of this class of objects. Four images of a distant (z s = 1.695) 
quasar are created by the gravitational field of a fairly nearby (zd = 0.0394) galaxy that is 
lying nearby in the line of sight to the observer. As a result, the time delay between the light 
variations for different images is less than a day, and the characteristic duration of the HMEs 
should be several tens of days. Uncorrelated fluctuations in the fluxes from the different images 
that were probably associated with microlensing by stars in the lens galaxy were first detected 
by Irwin et al. [9] and used to estimate the size of the accretion disk [10,11]. Later, various 
groups of observers monitored this object in hopes of detecting the effects of microlinsing [12-15]. 

The most complete series of observations, obtained by the international OGLE group [16,17], 
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demonstrated the presence of possible HMEs in components A (in 1998) and C (in Summer 
1999). The observational data for these events were analyzed with the aim of deriving the 
size of the accretion disk, both using statistical methods [18,19] and via model- fitting of the 
brightness distribution for a circularly symmetric source [20,21]. 

Alcalde et al. [22] recently presented observations by the GLITP group that nicely supplement 
the data of the OGLE group in Autumn 1999; these display the HME in component A. Based 
on the hypothesis that this event was associated with microlensing by a fold caustic, these ob- 
servations were analyzed by fitting both symmetric source models [23] and a model brightness 
distribution of the form expected for a standard geometrically thin and optically thick New- 
tonian accretion disk [24] surrounding a Schwarzschild black hole, allowing for the inclination 
of the plane of the disk to the line of sight [25]. In the latter case, it was possible to derive 
constraints on the mass of the black hole at the center of the quasar: 1O 7 M < M < 6xlO 8 M . 

Together with many advantages, the model-fitting also has certain disadvantages. The main 
one is the possibility that model is inadequate to the real object. It is clear, for example, that 
the possibility that plane of the disk is inclined to the line of sight makes circularly symmetrical 
models potentially inadequate. Similarly, there is a good basis to suppose that there should be 
a Kerr black hole rotating at close to the maximal rotation rate at the center of a quasar [1]. 
In this case, the accretion disk is appreciable relativistic. The edge of the rotating disk that is 
approaching the observer will appear to be brighter than the receding edge, leading to a loss 
of symmetry in both radial brightness distribution and strip brightness distribution across the 
disk. Therefore, model-independent methods for the analyses of observational data based on 
the restoration of the brightness distribution are especially important. 

The important censorious remark must be made relative to the use of the Tikhonov's regulari- 
zation method. Unlike the binary-lens case, in which caustics are isolated and their crossing 
by a source can be clearly fixed, microlensing by stars of a lens galaxy produces a random 
field of caustics in the source plane. The neighboring caustics also contribute to observed flux 
variations. Their influence can be minimized by analyzing only part of the light curve near 
its maximum. In this case, however, it becomes impossible to determine the initial flux level, 
relative to which the flux varies as additional images appear or disappear during the primary 
caustic crossing that led to the HME. This initial flux level cannot be derived directly from the 
light curve. 

At first glance, it seems that this problem has a simple solutions - allowing the initial level to be 
a free parameter and estimating it by requiring the best fit to the observational data. However, 
this approach automatically makes it impossible to apply the Tikhonov's regularization method 
used in [6-8] to restoration of the brightness distributions. It is known that the approximate 
solution of the ill-posed inverse problem yielded by this method is the smoothest function whose 
lensing curve agrees with the observational data within the errors [26,27]. Varying the initial 
level produces a family of strip brightness distribution. All functions of this family, including 
those having no physical meaning, will give the same goodness-of-fit to the observations. There 
is no criterion that can be used to select one of these function as the best approximate solu- 
tion. In addition, the Tikhonov's regularization method makes the minimum use of a priori 
information about the solution of the ill-posed problem, and is therefore not able to provide 
good stability of the solution at large noise level. 
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The purpose of our paper is the description of the new technique of the analysis of observational 
data by invoking additional a priori information consistent with the physics of the phenomenon, 
and results of its application to the analysis of the HMEs observations in components of the 
gravitational lens QSO 2237+0305. These results were published partially in the papers [28,29]. 

2. Method for restoration of strip brightness distribution across source 



Let us suppose that the image of a quasar's accretion disk is scanned by a fold caustic, which 
can be taken to have the form of a straight line due to the small angular size of the disk. Let 
b(x, y) be the brightness distribution in the disk for a distant external observer in a Cartesian 
coordinate system (x, y) in the plane of the sky, with the x axis oriented perpendicular to the 
caustic and the coordinate origin coincident with the center of the disk. The observed lensing 
curve will then depend only on the one-dimensional strip brightness distribution B(x), defined 
by expression 

/oo roo 
/ b(S,y)6(Z-x)dSdy, (1) 
-oo J —oo 

where 8{x) is the Dirac delta function. When the caustic crossing leads to the appearance of 
additional images of the source accompanied by a sharp amplification of the flux, the HME 
light curve I(x) is given by the convolution integral equation 

/oo 
A(x-Z)B(Z)d£, (2) 
-oo 

whose kernel has the form [30-32] 

A(x -t)=Ao + -J=l H{ * x - 0» (3) 

where H(x — £) is the Heavyside step function, which is equal to zero for negative and unity 
for non- negative values of its argument. The quantity A in (3) depends on the initial flux 
level, which is produced by all the other images of the source and remains unchanged during 
the caustic crossing, and the factor K characterizes the amplification of the caustic. The values 
of A Q and K are usually not known. When (3) is substituted into (2), we can see that A 
determines the initial flux level 1$ , which as we said in the introduction, is a free parameter 
of the inverse problem. The absence of information about K means that B(x) can be restored 
only up to a constant factor. The second free parameter is the time when the caustic passes 
through the center of the disk t , which determines the origin for the x axis. If the projection 
of the tangential velocity of the caustic onto this axis is V±, the time dependence of the spatial 
variable x has the form x — V±(t — to). As for the parameter K from (3), the scanning velocity 
Vj_ is, in general, unknown. Setting K = 1 and V± — 1 and assuming x = t — t , we can restore 
only the form of the strip brightness distribution B(x) from observations of the HME. 

Since the quasar accretion disk is appreciably relativistic, if the rotational axis is inclined to the 
line of sight, the brightness of an area element approaching the external observer should exceed 
the brightness of a receding area element [7]. Therefore, apart from cases where the caustic 
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scan is directed along the rotational axis, B(x) ceases to by an even function. We can now 
write the initial integral in (2) as a sum of two integrals taken along negative and non-negative 
intervals. Performing a substitution of variables in the first of these integrals, we obtain 



It can easily be seen that the properties of the kernel (3) imply that for x G (— oo,0] the 
second integral in (4) vanishes. It follows that only the negative branch of the strip brightness 
distribution contributes to the formation of the negative branch of the lensing curve, whereas 
both branches affect the positive branch I(x) when x G (0, oo). 

The analysis of the strip brightness distributions for relativistic accretion disks [28] shows that 
for optical and IR radiation outside the small region | £ |< £ > which makes a negligibly small 
contribution to the total flux, £?(—£) for £ G (—oo, — £o] and for £ G [£o, oo) are either non- 
negative, monotonically non-increasing or non-negative, monotonically non-increasing, convex 
downward functions. It is known that the sets of functions of these types are the compact sets. 
The search for the solution of an ill-posed problem on the compact set of functions gives the 
unique and stable result [26,27]. Thus, we can regard the branches of -B(£) as members of the 
compact set of these functions. This a priori information is qualitative and imposes no rigid 
model constraints on the form of the strip brightness distribution. 

The values of free parameters Iq and to can be determined from the minimum residual corres- 
ponding to the restored function B(x), which is usually found using the sum of the squared 
deviations, or the quantity Xn ■ Since the set of functions on which the search for the solution 
is carried out is compact, this guarantees that the obtained profiles B(x) and the values of the 
free parameters will approach their exact values as the errors in the estimate of the observed 
flux I(x) approach zero [26,27]. The use of a large amount of a priori information about the 
possible form of the strip brightness distribution in accordance with the physics of phenomenon 
enables us to achieve a solution with a high degree of stability against the effects of random 
noise. 

Given the above considerations, we have formulated the following algorithm for successive 
restoration of the branches of the strip brightness distribution assuming K — 1, V± — 1 and 
x = t — t . 

Step 1. Specify the initial values of the free parameters Iq and to • 

Step 2. On the negative branch of the lensing curve I(x) for x G (— oo,0], solve the inverse 
problem for the integral equation 



and find the negative branch B(— £) for £ G (— oo, — £o] 011 one of the compact sets of functions. 
Step 3. Modify the positive branch of the lensing curve for x G (0, oo) as follows: 



poo POO 

I(x) = A(x + C)B(-0 di + A(x- 0B(0 dC 



(4) 




POO 

I(x)= / A(x + 0B(-0dt 



Jo 




Step 4- Solve the inverse problem for the integral equation 
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for x E (0, oo) and find the positive branch of B(£) for £ G [£o, oo) on the same compact set of 
functions. 

Step 5. Compute the value of the total residual function for both branches of the lensing curve. 

Repeating the steps 1 - 5, we can easily find the global minimum of the residual function by 
exhausting all possible free parameter combinations. The values of parameters Iq and t and 
the two branches of the strip brightness distribution corresponding to this minimum yield the 
optimal approximate solution of the problem. 

A similar technique for the successive restoration of the branches of B(£) can be used when 
a caustic crossing involves the disappearance of source images, resulting in an abrupt flux 
decrease. In this case, the argument of the kernel (3) has the opposite sign. Therefore, only 
the positive branch of B(£) contributes to positive branch I(x), whereas both branches of the 
strip brightness distribution contribute to the formation of the negative branch I(x). 

We should note one other important circumstance, which has usually been neglected in the 
studies of the HMEs cited above. The main integral equation (2) has the singular kernel (3). 
When calculating such integrals, it is necessary to take special measures to ensure convergence 
of the corresponding integral sums. In particular, attempts to calculate the values I(xi) on non- 
uniform grid Xi when -B(^) is specified on a uniform grid can lead to large errors in the results. 
General questions with regard to the application of numerical methods for singular integral 
equations are considered in [33]. A simple proof of a sufficient condition for the convergence 
of the integral sums for the specific case of (2) with the kernel (3) is presented in [34]. This 
condition consist of the special selection of grids forming a so-called canonical dissection of 
the integration interval. Both grids are uniform, and either ^ = (xj + Xj + i)/2 or an integer 
number of grid steps A£ fit into the interval Ax, with the points £j being the centers of these 
intervals. Time span between measurements in the observed HMEs light curves are usually non- 
uniform. Therefore, the model light curves that are to be fit to the observations must first be 
calculated on a uniform grid, then interpolated to the moments of observations. When restore 
the brightness distribution, observed values I(xi) must first be interpolated to the uniform grid. 

3. Result of numerical simulations 



It follows that our technique for restoration of the strip brightness distribution is model inde- 
pendent. However, it is reasonable to test the potential of this technique by applying it to a 
realistic model of a quasar accretion disk. The theory of disk accretion onto compact objects 
initially developed by Shakura and Sunyaev [24] was further refined to include relativistic effects 
by Novikov and Thorne [35], and Page and Thorne [36]. We have used this standard model of 
an accretion disk. 

We considered the geometrically thin, optically thick accretion disk rotating in the prograde 
direction in the equatorial plane of a Kerr black hole with mass M = 1O 8 M , maximal 
normalized angular momentum a = 0.998 [37] and luminosity close to the Eddington limit. 
The details of these calculations are given in the paper [28]. The solid line in Fig.l shows the 
strip brightness distribution of our disk model B(x) in the V photometric band for the angle 
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between the rotational axis of the disk and the line of sight % — 45. The rotation of the disk 
is counter-clockwise and x axis is coincided with the major axis of the ellipse which is the 
projection of the disk onto the picture plane. The coordinate origin coincident with the center 
of the disk and the scale along the x axis is in units of the gravitational radius r g = GM/c 2 , 
where G is the gravitation constant and c is the velocity of light. We set the total flux emitted 
by the visible surface of the disk equal to unity. The negative branch B(x) in Fig.l, which 
correspond to the approaching edge of the disk, is clearly brighter than the receding side. Both 
branches - positive B(x) and negative as the function B(—x) can be approximated well by 
convex downward, monotonically non-increasing functions. 

In our numerical simulations, we used for calculation the HME light curve the first version 
of canonical dissection, with £j = (x, + x i+1 )/2, A£ = Ax = 20, & = —990 + 20(i — 1) for 
i = 1,2, ...M = 100 and x { = -1000 + 20(i - 1) for i = 1,2, ...N = 101 (both f and x in unit 
r g ). It be assumed that caustic crossing of the disk with B(£), shown by solid line in Fig.l, is 
accompanied by a flux decrease. We added the noise in the form of a random Gaussian numbers 
with zero mean and a standard deviation equal to 1% of the maximal flux value to the initial 
data samples I(xi). The resulting samples for the light curve I (xi) for the values of the free 
parameters Iq = and t = are shown by circles in Fig. 2. 

We searched for the branches B(£) and B(—£) on the compact set of non- negative, monotoni- 
cally non- increasing, convex downward functions for various Iq and t values using a modified 
version of the PTISR code written in FORTRAN [26,27]. To reduce the effect of roundoff errors, 
we transformed all real variables used in the PTISR and its auxiliary subroutines into double 
precision variables with 16 significant digits in their floating-point mantissas. We also used the 
additional a priori information that B(£) is equal to zero at the ends of the domain of variation 
of the argument. This is equivalent to specifying the size of this interval, which required in 
the solution of any ill-posed problem. Estimating the size of the domain where B(£) takes on 
nonzero values does not present difficulties in practice. Initially, a domain clearly exceeding 
probable values for this interval is adopted, and its size is then refined via a series of successive 
approximations. 

It is convenient to adopt the quantity 

N 

X 2 N = 7 °(^) ~ Ic{Xi))/0-iY, 
i=l 

for the residual, where N is the number of data samples on the HME light curve I (xi), Oi are 
the estimated standard errors of these samples, and I c {xj) are the fluxes corresponding to the 
restored strip brightness distribution. The residual had its minimum for the free parameters 
values Iq = 0.000 ± 0.001 and t Q = 3.0 ± 2.0 . The samples for the corresponding restored 
branches of the strip brightness distribution are shown by circles in Fig.l, and the light curve 
is shown by solid line in Fig.2. The minimum x% = 50.2 , whereas the value corresponding to 
a probability of 50% that the hypothesis in question should be adopted is 100.3 for N = 101 
degrees of freedom. On the whole, the inferred free parameter values end strip brightness 
distribution branches are close to their initial specified values. Thus, the results of our numerical 
simulation show that the proposed technique is potentially a powerful tool. 
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4. Observational data for HMEs in components of 
gravitational lens QSO 2237+0305 



The observational data of the OGLE group show that, in Summer of 1999, component C 
exhibited a characteristic flare, such as may occur during microlensing by a fold caustic accom- 
panied by the disappearance of additional source images. We obtained the V magnitudes for 
the HME light curve and their errors via the Internet from the OGLE server and transformed 
these magnitudes into flux samples averaging them over each observing night and assuming 
that a unit level corresponds to 18 m .O. These samples shown by circles in Fig. 3 as a function 
of time expressed in modified Julian date JD - 2450000.0. The vertical line segments indicate 
intervals corresponding to two standard deviations (±er). 

The HME in component A, observed in Autumn 1999, correspond to case when the caustic 
crossing leads to the appearance of additional images of the source accompanied by a sharp 
amplification of the flux. We obtained the V - band fluxes of component A (in millijansky, 
mJy) measured by the OGLE and GLITP (PSFphotll photometry) groups and their errors 
also via the Internet. The GLITP data, Ig(U), form a fairly dense series, but cover only the 
upper part oh the ascending branch and the maximum of the HME light curve. The OGLE 
data, Io(U), were more sparse, but cover the lower part of the ascending branch of the curve 
well. Therefore, we decided to use both of these photometric series in our subsequent analysis. 

When the OGLE and GLITP photometric data plotted together, it is obvious that there is an 
offset between them, and that the Io{ti) values exceed the Ig{U) values measured at the same 
moments. Since this offset is not large, and the data obtained in each photometric system cover 
a fairly extensive time interval with appreciable flux variations, we decided to unify the data 
series by assuming a linear relation between them: Io{U) = o>Ig(U) + b . The coefficients in 
this relation were determined from a least-squares fit with a linear interpolation of the Ig{U) 
values in the overlapping time interval, and were found to be a = 1.125 and b = —0.042 . 
Figure 4 shows the flux measurements reduced to the OGLE photometric systems as a function 
of the modified Julian date, JD - 2450000.0 . The filled and open circles show the OGLE and 
GLITP observations, respectively. The vertical line segments indicate intervals corresponding 
to two standard deviations (±er). The flux values do not display systematic differences, and 
the unified series of observations is fairly uniform. 



5. Results of restoration of strip brightness distributions 



Since the singular nature of the main integral equation (2) requires the use of grids with canoni- 
cal dissections, we again adopted the first version of such grids described above with Ax = A£. 
We searched for the branches B(£) and B(—£) on the compact set of non-negative, mono- 
tonically non-increasing, convex downward functions for various I and t values. Adopting 
K — 1, V± — 1 and x = t — to, we measure the distance x and the variable of integration £ in 
units of time. 
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The samples I (xi) were computed for the observations of the HME in component C via a 
linear interpolation onto a uniform grid x of N = 38 values with a step Ax = A£ = 6 d .O in the 
interval [to — 96 a! ,to + 126 d ]. The search for the branches of the strip brightness distribution 
was carried out on the grid & = — 237 d .O + (i — 1)A£, i = 1,2, ...M = 80. The values of the 
free parameter that yielded the global minimum of the residual were J = 1.020 ± 0.005 and 
t = 1392 d .3 ± d .l . The minimum of the residual function is Xn = 22.6, whereas the value 
corresponding to a probability of 50% that the hypothesis in question should be adopted is 37.3 
for N = 38 degrees of freedom. Figure 5 shows the restored branches of the strip brightness 
distribution. The corresponding light curve is shown by the solid line in Fig.3, and fits the 
observational data well. 

The canonical dissections with Ax = A£ = 2 d .O was used for the analysis of the observations of 
the HME in component A. The light curve was interpolated onto a uniform grid Xi of N = 69 
points in the interval [to — Q8 d ,to + 68 d ]. The search for the strip brightness distribution was 
carried out on a grid of M = 68 points in the interval [— 67 d , 67 d ]. The best-fit values of the 
free parameters were Iq = 0.697 ± 0.001 mJy and t = 1479 d .4 ± d .l. These are fairly close 
to the values obtained via model- fitting of the GLITP observations [23]. The values for the 
restored branches of the strip brightness distribution are shown by the circles in Fig. 6 and the 
light curve corresponding to the derived brightness distribution by the solid curve in Fig.4. 
The minimum residual was Xn = 37.2; the value corresponding to a 50% probability that the 
curve is agreement with the data for the case of N = 69 degrees of freedom is 68.3 . Thus, 
the agreement with the observations is fairly good. This is confirmed in Fig.4, where we can 
see that the calculated light curve tracks the observed counts well. At the same time, the 
curve is not completely consistent with the expected shape of the lensing curve expected for a 
fold caustic. This is especially noticeable in the interval following the flux maximum, where a 
deviation from monotonic behavior is observed. This was also pointed out in connection with 
model- fitting of the brightness distribution [23]. It is possible that curvature of the caustic, 
nearness to a cusp, or the influence of other nearby caustics is manifest in this HME. 

The presence of possible deviations from a simple model with a linear fold caustic is also 
reflected by the shape of the restored strip brightness distribution. We can see in Fig. 6 that 
two branches differ appreciable. The shape of the positive branch is close to the expected strip 
brightness distribution for a relativistic accretion disk ( see Fig.l). while the negative branch 
forms a simple straight line segment. Nevertheless, the derived strip brightness distribution is 
in agreement with the results of fitting the GLITP data with a standard model for Newtonian 
accretion disk, whose radius in time units proved to be 39 d .6 [23]. 

6. Estimations of sizes of accretion disk 



As was already noted above, if the time for the intersection of the source by the caustic is 
known, the linear size of the source can be derived using the information about projection of 
the tangential velocity of the caustic onto the axis normal to the caustic, V±. In addition, the 
time intersection of the accretion disk depends on the inclination of the plane of the disk to 
the line of sight and the direction of the motion of the caustic relative to the major axis of the 
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ellipse corresponding to the projection of the disk onto the picture plane. Therefore, even V± 
is known very accurately, we can determine only a lower limit for a the linear size of the disk 
from the strip brightness distribution. 

Let us project the spatial-velocity vectors of all the objects participating in the HME onto 
a plane perpendicular to the line of sight from the observer to the center of the accretion 
disk. We denote V s , V d) and V Q to be the two-dimensional projected velocities of the source, 
gravitational lens, and observer, respectively. As was shown in [5], the two-dimensional vector 
of the projection of the velocity of the caustic onto this plane V can be written as 

^ = V a Vd D s Vo Dds 

1 + z s 1 + z d D d 1 + z d D d ' U 

where z s and z d are the redshifts of the quasar and the gravitational lens, and D s , D d , and D ds 
are the angular diameter distances between the observer and quasar, the observer and lens, and 
the lens and quasar, respectively. For the parameters of the gravitational lens QSO 2237+0305 
(z s = 1.695, = 0.0394) and reasonable velocities for the motions involved, the second term 
in (5) is the determining one. Therefore, we can obtain for the velocity with which the source 
is scanned by the caustic 

V ± « , (6) 

l + z d D d ' yj 

where V± td is the projection of V d onto the local normal to caustic. 

The precise value of V± is not known. A statistical analysis of the peculiar velocities of galaxies 
yielded the mean value 663 km/s, while the range of the projected velocity V± id for 90% con- 
fidence interval was roughly 100 < V± id (km/s) < 1000 [25]. The distances D d and, especially, 
D s in (6) is depend on cosmological parameters. Taking into account possible variations of 
models of Universe from the usual (Qq = 1) to the model with a dominance of the vacuum 
energy consistent with modern observational data (fl = 0.3, A = 0.7), the interval of possible 
caustic-scanning velocities V± becomes 765 < V± (km/s) < 10548 [25]. 

If we formally adopt a value of the velocity in the middle of this interval, V± = 5600 km/s, 
then we can estimated the sizes of the accretion disk. The intersection time for the profile of 
the strip brightness distribution from the HME in component C (see Fig. 5) is roughly 300 d , 
which corresponds to a linear size of 1.5 x 10 cm, or 4.7 x 10~ 3 pc. This time from the HME 
in component A (see Fig.6) is 80 d , which correspond a size 3.9 x 10 15 cm, or 1.3 x 10~ 3 pc. The 
difference of the sizes is not very large. It is possible that this difference could be explained 
by either a difference of V±, or by the different direction of the caustic's motion relative to the 
major axis of the ellipse corresponding to the projection of the accretion disk onto the picture 
plane. 



7. Conclusions 



Our proposed technique for the successive restoration of the branches of the strip brightness 
distribution of a quasar accretion disk via analysis of observations of the HME makes it possible 
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to take into account a large amount of a priori information that consistent with the physics of 
the phenomena. This ensures that the resulting solution is accurate and stable against a random 
noise. Both branches of the strip brightness distribution can be derived from the light curve 
without applying any rigid model constraints. These properties are confirmed by numerical 
simulations and results obtained by applying the technique to real observational data. 

We have analyzed the HMEs in components C and A of the gravitational lens QSO 2237+0305 
observed by the OGLE and GLITP groups using our model-independent technique. This anal- 
ysis has yielded the estimations of the form of the strip brightness distribution across the 
accretion disk. The results for the HME in component C are consistent with the hypothesis 
that we have observed a scan of the source by a fold caustic. The form of the strip brightness 
distribution corresponds to our expectation for a relativistic accretion disk rotating around a 
supermassive black hole. In spite of a small value of residual x%, the features in the light curve 
for the HME in component A and the restored strip brightness distribution suggest that the 
disk was not scanned by simple linear fold caustic. It is possible that curvature of the caustic, 
nearness to a cusp, or the influence of other nearby caustics is manifest in the details of the 
light curve for this HME. 

The sizes of sources, measured in units of time, 300 d for the HME in component C and 80 d for 
the HME in component A, are consistent with the results of model-fitting. The crude estimates 
for the linear sizes are 1.5 x 10 16 cm and 3.9 x 10 15 cm, respectively. This difference could 
be due to a difference in either the scanning velocity, or the direction of the caustic's motion 
relative to the major axis of the elliptical projection of the quasar's accretion disk onto the 
picture plane. 
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Figure 1: Branches of the strip brightness distribution for the model of a relativistic quasar's 
accretion disk (solid line) and the strip brightness distribution samples restored in the numerical 
simulations (circles). The scale along the x axis is in units of the black hole gravitational radius 
r g = GM/c 2 . 
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Figure 2: Samples for the model HME light curve with noise superimposed (circles) and the 
curve correspond to the restored strip brightness distribution (solid line). The units along x 
axis are the same as in Fig.l. 
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Figure 3: Samples of the HME light curve observed by the OGLE group in component C of 
the gravitational lens QSO 2237+0305 in V band (circles) and the light curve correspond to 
the restored strip brightness distribution (solid line). 
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Figure 4: Samples of the HME light curve observed by the OGLE (filed circles) and GLITP 
(open circles) groups in component A of the gravitational lens QSO 2237+0305 in V band and 
the light curve correspond to the restored strip brightness distribution (solid line). 
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Figure 5: Branches of the strip brightness distribution across the quasar's accretion disk restored 
from the analysis of observations of the HME in component C. The scale along the abscissa 
axis is in units of time (days). 
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Figure 6: Branches of the strip brightness distribution across the quasar's accretion disk restored 
from the analysis of observations of the HME in component A. The scale along the abscissa 
axis is in units of time (days). 
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